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The paper shows how the known, exact results for the two electron bound states can modify 
the ground state phase diagram of extended Hubbard model (EHM) for on-site attraction, intersite 
repulsion and arbitrary electron density. The main result is suppression of the superconducting 
state in favor of normal phase for small charge densities. 
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I. INTRODUCTION 

The Hubbard model appears in almost all areas of solid state physics. Its universality is connected with the 
fact that it describes both band movement of charges as well as local and nonlocal ~ in the extended model - 
correlations 0]. Treating its parameters as effective ones, the model has been used in research of magnetism, 
superconductivity and especially high temperature superconductivity (HTS), other various phenomena in the 
solid, charge orderings, phase separation etc., in materials like high temperature superconductors, bismuthates, 
Chevrel phases, amorphous semiconductors, heavy fermion materials, systems with alternating valence to name 
the few (for a review see, e.g., Ref. "^I). 

Unfortunately there are not many exact results concerning this model. Usually the results are obtained in 
specific limits: infinite dimensions, one dimension (the most numerous group), infinite repulsion, half-filled band 
or other specific band fillings. We can mention the exact solution in one dimension {d — 1) obtained by Bethe 
ansatz 3j, Lieb's ferrimagnetism Nagaoka ferromagnetism in repulsive, half-filled systems with one hole 
fiat band ferromagnetism , some bounds on correlation functions 0, S and a finding of Randeria ^3 1 
according to which, existence of two-electron bound states of s-wave symmetry is necessary and sufficient 
condition for appearance of s-wave superconductivity in d — 2 systems with low electron density. Let's also 
note that the mean-field BCS equations for superconductivity in the Hubbard model with effective attractive 
interaction between electro ns, in the limit of vanishing electron density turn into Schrodinger equations, which 
can also be solved exactly . 

A phase diagram in two dimensions for arbitrary n, a case of special interest due to the possible connection 
with high temperature superconductivity, is still under examination. The results obtained in the mean-field 
approximation show competition of phases: in half-filled band charge density waves (CDW) for > 0, super- 
conductivity for [/ < and spin density waves for f/ > P|. The calculations for n ^ 1 suggest possibility of 
phase separation for [/ < 0, 0| : electron droplets for large enough < and phase separation of CDW 
(with n = 1) with SS for > (PS[CDW/SS]), for n around half-filled band competing with the pure SS state 
in low density limit. 

This paper shows, how the known solutions for the bound states (including exact solutions of Schrodinger 
equation) can be used for modification of the ground state phase diagram of the extended Hubbard model. 

II. HAMILTONIAN AND THE SUPERCONDUCTING STATE 

We begin with the extended Hubbard Hamiltonian in standard notation: 

ija i ij i 

where we sum over nearest-neighbors (nn) sites only. U and W are treated as effective parameters. We use 
broken-symmetry Hartree-Fock approach (for details see Ref. ill). As we are interested in the properties of the 
superconducting state, we introduce averages of operators C-kiCk] in Wick's- type decoupling J. 4] of the four- 
operator terms in the Hamiltonian. Non-zero average of such pair-creating operators means phase-coherence 
among pairs, i.e. a presence of superconducting state, and serves as an order parameter (see Eq. (|2l)- 

Ho = Y^i^k - -p)cl,Cka + X(^fe4TC-fei + h.c.) + C , (2) 
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where: 



and /I = M — + zW)n, where z is coordination number of hypercubic lattice, Wk — W^k, £k = —t"fk, 
7;. = 2^^cosfcQ., a e {x,y,z). After diagonahzation of the Hamiltonian Eq. (jSJ we obtain quasiparticle 
energy: 



E, = ^{e,-fir + \A,\^ . (4) 

and a self-consistent equation for the gap: 

^k = ^Y.iWk-, + U)^te.nh^, (5) 

where /3 = l/fc^T, T is temperature and ks Boltzmann constant. The constant C in the Hamiltonian Eq. ^ 
can be expressed now as: 

The pairing potential in the singlet channel (see Eq. Q) takes on the separable form for the square lattice and 
nn interaction: U + Wki-k2 = U + W^kilk2/ z (retaining only the terms of s-wave symmetry), and that makes 
possible solving Eq. (j^J by an ansatz: 

Afe = Ao + A.^7fc , (7) 
what leads us to the set of self-consistent equations: 

A7 = - — ^E ^9(^0 + ^9^7)^9 > (9) 
n-l = -^Y.^e,-T1)F,. (10) 

where Fq = (tanh ^^^)/2£'g. In the case of the rectangular density of states (DOS) and pure on-site pairing 
we can obtain analytical solutions [l^l; in the case of the extended s-wave superconductivity (Eqs © - Cni) 
analytical solutions exist in the limit of low electron density. Introducing a new parameter: A-y/Ap, we can 
expand Eqs ||HJ) - pUI) . treating Ao as a small parameter. As a result we obtain a formula for critical value for 
appearance of superconductivity for given U and n in the ground state : 

^- = M(l-n)+8^/~2MV^ ' "^^^^ ' = 

and D = zt \s half-bandwidth unit. In the case of rectangular DOS this formula reduces to [T^ : 

Wcr{l + {n-lf{l + m/U))^ -At . (12) 

We can go a step further in our mean- field analysis and include Fock term p = jjYlika lk{c],^Cka) into calcu- 
lations. In Eqs (O, Q, ijHl - Cnj, £k must be changed into Sk = £k{^ + pW/zt) then, and we have to solve 
Eqs © - 11U|1 self consistently with the equation for the Fock term: p = J^k ^klkFk- Equations - (|12|l 
remain valid, with the change X X/{1 +pW/zt) where X — n, Wcr and U. 



III. LOW DENSITY LIMIT 



Going back to Hamiltonian Eq. Q we can obtain exact results in the low density limit. In the center-of-mass 
coordinate system we can expand the wave function of the two-electron bound pair in the basis of plane 
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waves (i.e., eigenstates of the hopping part of the Hamiltonian Eq. (^). We can easily find the equations for 
the coefRcients of the expansion, what finally yields a set of self-consistent equations for the wave function in 
the position space, in terms of lattice Green functions PHvj: 

^(r) = J2 G{E, P, r, r')g{r')i:{r') , (13) 
r' 

where G is lattice Green function defined by: 

G(E.P,r.r')^-Y.^^^. (14) 

and g{r) is diagonal interaction matrix, consisting of elements U and W. Eigenenergy equation takes the form: 

det[g - g-'] = , (15) 

where G is a matrix with elements Qij = G{E,P ,ri,rj). This is an analogue of Eqs (O - Let's note that 
in the case of the two-electron bound pairs A = and the role of the binding energy is played by 7^/2. For the 
hypercubic lattices these equations were solved and it has been found out that in one and two dimensions pairs 
for = bind for any negative U, while in three dimensions there is critical value for W J^. The formula 
for Wcr in the case of two-electron bound state reads : 



\Wcr 

2t 



2D 



+ ic~ir\ (16) 



where C — 1/-/V^j,(l — Jk/z)^^ is the Watson integral. This is an exact result. Remembering that C is 
divergent for lattices of dimensions d — 1 and d = 2 we can see, that for n — » /Lt/Z? — > —1, / — > (Eq. 
and Eq. (|16|l is a limiting value of the formula Eqs 1)11(1 and l|12l) for d — 1 and d = 2, as it should be. Not for 
d = 3 though; this case will be discussed later on. 

Let's note that Eqs ((ll|l . ((12|l and ((16(1 . are valid for any combination of signs of U and W and for large 
enough U < and < there is a second branch of solutions 16, 19] . The two branches are the two solutions 
which realize in the two opposite limits: U — +oo and U — — oo (or W = ±oo). The formal equations and their 
solutions in both these limits are the same, despite completely different physical situation. Nevertheless these 
are the specific cases of two distinct solutions. 



IV. RESULTS AND DISCUSSION 



In view of the Randeria's notice, described in the Introduction, in the case of s-wave symmetry in two 
dimensions we can use condition for existence of bound states as a condition for existence of superconductivity. 
In Fig. 1 the boundaries expressed by Eq. 1(11(1 for different lattice dimensionalities and electron densities are 
shown. For parameters U, W belonging to the area above the plotted lines (mostly in the 1st quarter of 
coordinate system) two-electron bound states, and what follows s-wave superconductivity in two dimensions, 
can not exist. The curves for n = in all dimensions are exact results. 



1.0 



-0.5 



-1.0 

















\ \\v 

\ \ \\ \ 

• \ ■■□ A o 

__s 












\ 








1 1 1 1 1 1 1 






-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
U/D 



FIG. 1: Critical values for existence of bound pairs and 
superconductivity for: n = for d — 3 (black circles 

- bound pairs only), n = and rectangular DOS (full 
line), n — 0.25, rectangular DOS and Fock term (dotted 
line), n = 0.25 and rectangular DOS (squares), n = 0.25 
for d — 2 (triangles) and n — 0.25 for d = 3 (white circles 

- superconductivity only). Line n — for d = 2 is the 
same as for the rectangular DOS. Half-bandwidth unit 
D — 4t for rectangular DOS and for d = 2 while D = 6t 
for d = 3. 



As n gets larger, the area of existence of bound pairs increases for W > and U < and decreases in the 
part of the diagram with W < and U > 0. All curves except the one for n = in d = 3 go through the point 
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with coordinates (0,0). This illustrates the fact that ior W = infinitesimally small U creates bound state 
in d = 2, while the threshold exists in d = 3. Nevertheless we do not have threshold in d = 3 for n 7^ - in 
agreement with Randeria's notion about necessity of bound states for superconductivity only in d = 2. Let's 
note that for large U and W the curves approach the asymptotes - for curves crossing through the axes origin 
the asymptotes are given by the formulas: Uas/t = — 16(71 — 1)^/(1 + {n— 1)^) and Was/t — —4/(1 + (n — 1)^). 
This is connected with the fact that in the 3rd quarter of the coordination system for U < and < there 
exist second branches of the solutions. As they have higher energy than the solutions described in Fig. 1 they 
do not modify the ground state phase diagram and are not shown here. 

The dotted line in Fig. 1 (and in Fig. 2) describes the results of calculations with inclusion of the Fock term, 
using Eq. 1)12(1 modified by the (1 + pW/zt) term, as was described in the end of Section 2. To simplify the 
calculations the Fock term from the normal state was used: p = n{2 — n). For W > Q this term broadens the 
band moving the system more into weak coupling limit and enlarging the normal state area (opposite behavior 
for W <Q). The same effect can be seen in Fig. 2. 

In Fig. 2 the boundaries of existence of bound states, Eq. H12|l (black symbols), are plotted on a ground 
state phase diagram of the extended Hubbard model for U < and W > . fo r arbitrary n and rectangular 
DOS, together with the phase boundary PS[CDW/SS]/SS taken from Ref. [l3 (white symbols). Above the 
lines with black symbols s-wave superconductivity can not exist in d = 2. This way the superconducting state 
is suppressed and normal state (NO) area is introduced into the phase diagram. Let's note that also the phase 
separated state PS[CDW/SS] is "reduced" to the NO phase and not to the CDW phase. This is due to the 
fact that the CDW in the PS state is the CDW with n = 1. CDW with n 7^ 1 is unstable, as it has negative 
compressibility. 
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FIG. 2; Phase boundaries for U/D = —0.4 (circles) 
and U/D — —0.8 (squares) calculated for the rectan- 
gular DOS. Black symbols denote boundary of existence 
of bound state, white symbols boundary between sin- 
glet superconductivity (SS) and phase separated area 
PS[CDW/SS] (with charge density wave CDW and SS). 
Symbols on dotted lines show the results of calculations 
including Fock term. 



Another thing to note is the threshold for appearance of bound states for n = 0, which increases with 
increasing \U\. The phase diagram is modified only for intermediate values of |?7| and \W\, smaller from their 
asymptotic values \Uas\ and |VFqs|. For \U\ or \W\ larger than these values, bound states exist for arbitrary 
value of the other parameter, in agreement with Fig. 1. 

The calculations in Ref. jl^l consider only pure, on-site s-wave pairing. Including A-y (Eq. 0) into calculations 
does not change much the described PS[CDW/SS]/SS boundary - A-y is two orders of magnitude smaller than 
Aq on this boundary. Inclusion of Fock term into the calculations of the bound states, results in extending the 
area of the normal phase, as was mentioned before. This effect increases with increasing \n\ and W. 

In conclusion it was shown, how the analytical (and exact for n = 0) formulas for bound two-electron states 
can be used for the modification of the U < 0, > part of the phase diagram of the extended Hubbard 
model. The main result of this approach is suppression of the superconducting and phase separated areas in 
favor of the normal phase around half- filled band for intermediate values of |C/| and \W\, larger than threshold 
values and smaller than \Uas\ and |W^as|- 
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